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ABSTRACT 

Recent mid-infrared observations of young stellar objects have found significant variations possibly 
indicative of changes in the structure of the circumstellar disk. Previous models of this variability 
have been restricted to axisymmetric perturbations in the disk. We consider simple models of a non- 
axisymmetric variation in the inner disk, such as a warp or a spiral wave. We find that the precession 
of these non-axisymmetric structures produce negligible fiux variations but a change in the height 
of these structures can lead to significant changes in the mid-infrared flux. Applying these models 
to observations of the young stellar object LRLL 31 suggests that the observed variability could be 
explained by a warped inner disk with variable scale height. This suggests that some of the variability 
observed in young stellar objects could be explained by non-axisymmetric disturbances in the inner 
disk and this variability would be easily observable in future studies. 
Subject headings: stars:pre-main sequence, circumstellar material, star:variables:other 



1. INTRODUCTION 

Variability has been one of the prominent features of 
T Tauri stars since they were first observed (Joy 1945J). 
In the optical fiux changes up to 0.5 magnitudes over a 
few days are common with most of the variability ex- 
plained by hot and cold spots r otating across the surface 
of the star (jHerbst et al.lll994 ). In a few cases, circum- 
stellar disk variability has been inferred from stochastic 
or quasi-periodic optical occultation events. In the case 
of AA Tau a warped inner disk occults the star lead - 
ing to periodic optical variations (jBouvier et al.ll2Q03l ). 
For HH 30, putative structural changes in the disk may 
lead to the observed changes in th e scattered light emis- 
sion in as little as a few days (jWatson fc Stapelfeldtl 
|2007| ). Variations in the scattered light images of HD 
163296 are also attributed to a cha nge in the structure 
of the disk (jWisniewski et al.ll2008[ ). FU Ori outbursts, 
which show an orders of magnitude increase in the op- 
tical flux followed by a slow decrease over decades, are 
associated with changes in the accretion rate of the disk 
brought about by thermal or gravitational instabilities 
(jHartmann fc Kenvon|[T996l ). EXor stars also show vari- 
ability on month to yearly timescales probably rel ated to 
changes in the accretion rate through the disk (jHerbiel 

Hop). 

Infrared observations provide a more direct way of 
looking for variability of the circumstellar disk. Emis- 
sion at these wavelengths comes from dust near the 
surface of the central <10 AU of the disk and is sen- 
sitive to changes in the disk temperature and struc- 
ture. The dust temperature is set by a combination 
of viscous energy release as disk gas accretes toward 
the central star and irradiation of the disk surface by 
light from the stellar photosphere and accretion shock. 
Variations in the near-infrared cannot always be ex- 
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plained simply by a change in the stellar flux, but a 
change in the st ructure of the c ircumstellar disk must 
also be included (ICarpenter et al|2 001: Eiroa et al. 200^ 
iSkrutskic et al."1996' 'A lves de O'liveira fc Casalil 12008^ . 
Plavchan et al. (2008) flnd variable disk emission at 3- 
8 /im around the star WL 4, which also exhibits pe- 
riodic near-infrared changes suggestive of influence by 
a binary companio n. Variabihty a s far out as lOO^tm 
has been observed (iLiu et al.lll996t iJuhasz et al.l 120071 : 
iBarsonv. Ressler. fc Marsh! 120051 ) although the limited 
wavelength coverage and time sampling do not provide 
strong constraints on the physical cause. Some of this 
variation can be explained by changes in the structure 
of the inner disk. Variability of the 10/im silicate fea- 
ture, which is sensitive to grain properties at the surface 
of the inner few AU of the disk, has also been observed 
suggesting that the grain properties of disks are rapidly 
varving ()Abraham et al.ll2009l: iSkemer et aLll2010D . 

Most disk models assume axisymmetry, and thus ex- 
plore axisymmetric structure variations or asymmet- 
ric irradiation fields to explain the observed variability. 
iMuzeroUe et al.l ()2009f) suggest that a non-axisymmetric 
inner disk may explain the strange wavelength depen- 
dence of the variations seen in LRLL 31, a transition 
disk in the young cluster IC 348. In this paper we further 
explore the possibility that these mid-infrared variations 
could be caused by a non-axisymmetric structure in the 
inner disk. First we consider three simple models for a 
non-axisymmetric structure: a warp in the middle of the 
disk, a warp at the inner edge of the disk and a spiral 
wave in the disk. We examine how the SED of these 
disks varies as the structures precesses/corotates around 
the star and as the height of the structure changes with 
time. We then apply these models to the variations in 
LRLL 31 as well as to other mid-infrared variable stars. 

2. MODELS 

2.1. Motivation 

A number of authors have recently found strong wave- 
length dependent mid-infrared variability in both Her- 



big AeBe stars and T Tauri stars. iJuhasz et al.l (|2007[ ) 
observe 3-100/im variability of SV Cep over a two year 
period. These variations can be up to a factor of two at 
lOO/im, and there appears to be a correlation between 
the lOO/im flux and the optical flux and a weak anti- 
correlation b etwee n the 3.6/im flux and the optical flux. 
iSitko et"aLl (|2008[ ) find significant variations in the 3- 
13/im flux of the Her big AeBe star HP 3 f 648 over almost 
a quarter-century. iHutchinson et al.l 0994) find two 
stars, UX Ori and AK Sco, whose optical fiux decreases 
as th e f 0/im fiux increases. Morales-Caldcron et al] 
(j2009( ) find that 3-8/im variability of 0.05-0.2 mag is com- 
mon in pre-main sequence star s in the young cluster IC 
1396A. iMuzeroUe et al.l ()2009[ ) find mid-infrared varia- 
tions of up to 60% in as little as one week in the young 
star LRLL 3f . The IRS spectra show a clear wavelength 
dependence of the fiux variations, where the 5-8//in fiux 
decreases and the 8-40/im fiux increases. This variation 
is especially interesting since LRLL 31 exhibits an SED 
typical of a transition disk, indicative of a deficit of small 
grains in the inner disk. 

Typical sources of optical variability in pre-main se- 
quence stars have difficulty explaining the observed mid- 
infrared variations in these stars. Cold spots on the 
star cannot produce a large enough variation in the mid- 
infrared, where the firrx is dominated by the circumstel- 
lar disk. They would not be able to produce variability 
out as far as 100/im, as has been seen in a number of 
stars. Hot spots, due to the accretion fiow striking the 
surface of the star, produce large optical variations that 
may affect the infrared fiux. As the star rotates, the hot 
spot acts as a fiashlight heating a localized portion of 
the disk. If the disk is viewed close to edge-on then the 
asymmetric heating can produce noticeable variations in 
the infrared (Morales-Caldcron ct al. 2009). However, it 
is difficult for this model to explain some of the larger 
variations and different wavelength dependences seen in 
many young st el lar ob jects. 

iJuhiisz et ail (|2007[ ) and ISitko et al.l (|2008D use de- 
tailed radiative transfer models of a varying height at 
the inner edge of the disk to reproduce their observations. 
Increasing the scale height raises the near-infrared fiux 
while shadowing the outer disk, and occulting the star 
if the system is close to edge on. The long-wavelength 
variation is difficult to explain with in-situ changes in the 
structure of the disk because the dynamical timescale 
is years to decades longer than the observations, thus 
models with variations in the inner disk that shadow the 
outer disk are used. Due to the restrictions of the ra- 
diative transfer code, the perturb ations were assu i ned to 
be axisymmetric within the disk. IMuzeroUe et al.l (|2009[ ) 
suggest that the variation in LRLL 31 may be a warp, 
spiral wave or some other non-axisymmetric structure in 
the inner disk. Some of the observations of LRLL 31 
are difficult to explain with a puffed inner rim and other 
configurations were considered. Here we follow up on 
that possibility and develop a simple model of such struc- 
tures. When examining these non-axisymmetric models 
we consider two time-dependent effects: azimuthal mo- 
tion such as precession and corotation, and variations in 
the scale height of the perturbation. Precession is a nat- 
ural consequence of the torques that originally warped 
the disk while a variable scale height may be related to 
the dynamical process that causes a warp in the first 



place. Both of these are possible sources of variability 
for a non-axisymmetric structure. 

2.2. Deriving the SED 

Disks are generally modeled as axisymmetric 
fiarcd disks that are heated by stellar irradiation 
(Chiang & Goldrcich 1993) or a combination of stellar 
irrad iation and viscous dissipation from accretion 
(D'Alessio et al.l I2OO60 . But the disks can become 
disturbed from this initial state by the gravitational 
influence of a binary com pa nion or passing star 
(lArtymowicz fc Lubowl 119941 : iLarwqod fc Papaloizq 
19971 or an instability within the disk (jWatanabe fc Li: 
2008) . This can lead to the disk becoming warped, which 
will change the temperature structure and emission of 
the disk. To study how structural changes will effect 
the SED we consider three simple models (1) a disk 
with a middle warp, (2) a disk with an inner warp and 
(3) a disk with a spiral wave, shown in figure [TJ All of 
these disks are taken to be geometrically thin blackbody 
emitters. Flaring and an optically thin atmosphere 
will have a substantial effect on the SED but these 
simple models present a useful starting point. If our 
models were completely confined to the midplane the 
SED would have the simple form \Fx ex \~'^^^. How 
our models depart from the power law will depend on 
the location of the structural disturbance within the 
disk. To calculate the SED for o ur simple models we 
follow the derivation presented in iTerguem fc Bertout] 
(11996', hereafter TB96). TB96 lay out a procedure for 
estimating the temperature structure of a blackbody 
disk whose height, h, is an arbitrary function of r and 
9 within the disk. The coordinate system is defined as 
centered on the star with 9 = along the line with the 
highest warp above the midplane. The fiux intercepting 
a point P{r, 9) in the disk is given by: 
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with tan{/3) = h{r,9)/r. The variables S and e are the 
polar and azimuthal coordinates of a point on the surface 
of the star with respect to the axis joining the polar axis 
to the center of the star. Converting to temperatures, 
and performing the integration over e, the equation then 
becomes: 



T\r,i 



— —T^ 



2.3. Middle Warp 
We first consider a warp similar to that originally dis- 
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The last term in this equation is the energy from vis- 
cous dissipation, which for typical T Tauri accretion rates 
is much less than the heating from irradiation. We as- 
sume that the accretion rate onto the star is small enough 
that the accretion luminosity is negligible compared to 
the stellar luminosity and thus does not contribute to 
the heating of the disk. The process of evaluating the 
temperature of the disk reduces to a determination of 
5max,Smin,£min,£max- The three modcls win vary in 
their evaluation of Smax,Smin,£min,£max for each point 
P{r, 9) on the disk surface, which is discussed in more 
detail in the appendix. 

From TB96, the flux emitted by the disk is: 
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In this case u is the vector along the line of sight to 
the observer from the center of the star and nd is the 
normal to the disk at the point P{r, 9). The vector u can 
be defined in terms of the azimuthal and polar angles to 
the line of sight, a and i respectively. The disk surface 
extends from rmin, which we take to be 8i?*, out to Rdisk, 
typically 100 AU. The inner radius 8i?* is roughly equal 
to the dust destruction radius for a frontally illuminated 
disk with wel l-mixed small dust gra ins around a 0.5Mq 
T Tauri star (jMuzeroUe et al.ll2003l ). 

The area of the disk along the line of sight is given by 
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We also include the possibility of the star being oc- 
culted by the disk. These prescriptions are described 
in detail in the appendix. This is highly dependent on 
the viewing angle, which for the non-axisymmetric disks 
considered here depends on both inclination, i, and the 
azimuthal viewing angle, a. The line of sight a = 
corresponds to looking at the disk along the line 6 = 0, 
while an observer along i = 0° is viewing the disk face- 
on. TB96 do not examine the dependence on azimuthal 
angle, although their generic equations repeated above 
do include a. In the appendix we demonstrate how this 
dependence has been included in the models. 

For our fiducial models we take the central source to be 
a G6 star with a radius of 2.9 Rq, a mass of 1.8 Mq and 
an accretion rate of M = 10~^MQyr^^. These specific 
parameters are chosen in anticipation of modeling the 
disk around LRLL 31. As we will show later, the shape 
and strength of the variability does not depend strongly 
on the central star and the results presented in the next 
sections are applicable to a broad range of pre-main se- 
quence stars. 
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for r ^ Rwarp- The parameter Rwarp specifies the 
outer radius of the warp. TB96 take Rwarp to be equal 
to the radius of the disk, while here we consider situa- 
tions where the disk extends beyond the warp, shown in 
figure la, in order to examine shadowing effects. Past 
Rwarp the disk is confined to the midplane and this flat 
passive extension can be shadowed by the warp, chang- 
ing its temperature structure. Since the warp occurs in 
the middle of the disk we refer to it as a middle warp. 
We consider disks with Rwarp of a few AU or less in order 
to probe any effects on the 2 to 24/im flux. The param- 
eter g sets the height of the warp, and is equivalent to 
the maximum value o f fe/r within the disk. According to 
iTerguem fc Berto"utl (|1993l ) a warp created by an equal 
mass binary will have n = 4, g = 0.01. Figure [5] shows 
the temperature structure of disks with g = 0.01,0.05 
and Rwarp = 0.5 AU. These warps demonstrate the shad- 
owing of the outer disk and how it can change with ra- 
dius. Immediately behind the warp (dark line beyond 0.5 
AU) the disk is completely shadowed and is only heated 
by viscous dissipation. Further out the warp only blocks 
part of the star and the temperature approaches that of a 
flat passive disk. The height of the warp will change the 
temperature at the warp as well as the shadowing of the 
outer disk. Figure [2] compares the SED of this warped 
disk with a flat passive disk and a flat accretion disk. 
The short wavelength flux rises above a flat passive disk 
because of the high temperature on the side of the warp 
facing the star. The shadowing of the outer disk reduces 
the flux at long wavelengths so that it is less than a fiat 
passive disk. 

A warped disk will be driven to precess by the torques 
on the system which originally warped the disk. To 
mimic a precessing warped disk, we vary the azimuthal 
viewing angle, a, of the observer. The method of TB96 
needs to be modified to account for a changing azimuthal 
angle, and these changes are described in the appendix. 
In figure [3] we show the precession of a warp with 
g = 0.05 which extends to Rwarp = 0.5AU with a fiat disk 
extending to 100 AU. As the disk precesses the infrared 
excess only changes by ~ 0.05 dex. The flux change at 
short and long wavelengths is correlated. The observed 
wavelength dependence can be understood by consider- 
ing the schematic middle warp drawn in figure [TJ The 
observer is rotating from the upper right to the upper left 
in this drawing and as the observer rotates the hot side 
of the warp comes directly into view. This increases the 
short wavelength flux, and since the orientation relative 
to the cold outer disk does not change with a the long 
wavelength flux does not change significantly. The dis- 
tance between the warp and the star keeps the warp cold 
and causes the change in the SED to be at longer wave- 
lengths. Figure [3] also shows the precession of a warp 
seen at « = 85°. The infrared excess is smaller, because 
of the more highly inclined observing angle, while the 
change in fiux is larger although the wavelength depen- 



dence is the same. As the disk becomes more inclined the 
warp is seen closer to face-on, creating more significant 
variations as it rotates in and out of view. 

Creating a warp in the middle of the disk is possi- 
ble with a (sub)stellar companion embedded in the disk. 
If the companion is not coplanar with the disk then the 
gravitational perturbation from the companion can drive 
material out of the midplane. The disk around /? Pic dis- 
plays a war p in the middle of the disk due to an unseen 
companion (jAugereau et al.l[2001l ) . The variable star WL 
4 displays periodic variability in the near-infrared, in- 
dicative of a multiple system, as well as mid- infrared vari- 
ability, i ndicative of a change in the dust structure near 
the star (jPlavchan et al.ll2008D . Roughly 10% of stellar 
syste ms have a binary companion within 1 AU (jMathieu I 
Il994f ). If the companion is within one AU of the central 
star then the precession period of such a configuration is 
years to decad es, depending on the exact configuration 
of the system (jCaprani et al.|[2Q06[ ). The combination 
of the small flux change and long timescale would make 
precession difficult to observe. 

We next consider the effect of changing the height of 
the warp. As the warp height increases the disk will get 
heated to a higher temperature since it is more directly 
illuminated by the radiation field. This is the same rea- 
son that a flared disk gets heated to a higher tempera- 
ture and has a higher flux than a flat passive disk. The 
larger warp will also cast a larger shadow on the outer 
disk which will decrease its temperature and lower its 
flux. If the warp changes substantially enough then the 
change in the temperature structure would be evident in 
the SED. The hot surface layers of the disk that are seen 
in the mid-infrared respond almost instantaneously to a 
chang e in the stellar radiation field ([Chiang fc GoldreichI 
Il997f ). This means that the timescale for the SED to 
change will depend on how quickly the disk can change 
its structure, which is the dynamical timescale. For the 
central star considered here the dynamical timescale can 
be as low as one week within 0.1 AU. 

In order to examine this effect we take the disk we 
used to study precession {Rwarp = 0.5AU) and vary its 
scale height, parameterized by the factor g, from 0.01- 
0.1. This corresponds to a physical height of 0.005-0.05 
AU. The results are shown in figure U) We find that as 
the height of the warp increases the flux at A < 15/im in- 
creases while the flux longward of 15/im decreases. The 
flux variations is small but the unique wavelength depen- 
dence would make infrared variability due to a changing 
scale height at the inner edge of the disk easy to diagnose. 

The same mechanism that drove the disk to become 
warped may also cause it to change the height of the 
warp. If the warp is caused by a companion, then this 
object cou ld drag material as its or bit takes it out of the 
midplane (jFragner fc Nelson! I2009f ) . For a binary com- 
panion whose orbit is misaligned with the disk, the plane 
of the disk will equilibrate with the plane of orbit over 
thousands of years, but on the timescale of a single or- 
bit dust will be shifted by the passing of the companion. 
As the companion's orbit takes it out of the plane of the 
disk it will drag material with it, but when the compan- 
ion's orbit takes it back into the midplane of the disk this 
dust will quickly settle down. This would lead to peri- 
odic variations in the scale height of the warp on top of 
the long term variations. It may also be possible that a 



long-lived asymmetric structure within the disk, such as 
a warp caused by a companion, is experiencing variable 
illumination from the central source. This is possible if 
the accretion flow is localized into hot spots that rotate 
around the star. As the hot spot illuminates the warp, it 
will heat the side facing the star, causing its scale height 
to increase. The scale height of the warp will decrease as 
the hot spot rotates around to the far side of the star. 

2.4. Inner Warp 

As the warp moves closer to the star, the deviation in 
the SED from a flat passive disk moves to shorter wave- 
lengths. To affect the flux at 5— 8/im we need to move the 
warp almost to the inner edge of the disk. We could also 
consider a warp that reaches its maximum height above 
the midplane at the inner, rather than the outer, edge 
of the disk. An inner warp has bee n invoked to explain 
the variability observed in AA Tau ([Bouvier et al.|[2003l) 
in which the stellar magnetic field is misaligned with the 
disk causing it to warp as material flows onto the fleld 
lines. A inner warp can be described by the function: 
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The temperature distribution for disks with g — 
0.05,0.01, rjnin = 8R* are shown in flgure [5j The ex- 
ponent for the power law was chosen to match the mid- 
dle warp and depends on the exact physical cause of the 
warp. The inner warp can be split into two pieces, a 
convex piece which directly faces the star and a con- 
cave piece which is turned away from the star. In the 
schematic drawing of the inner warp in figure [T}d on the 
right side the convex piece is on the bottom while the 
concave piece is on the top. The inner edge of the con- 
vex side reaches a very high temperature and departs sig- 
nificantly from a flat passive disk. The concave side just 
behind the warp has a low temperature because the warp 
completely blocks the star and the disk is only heated by 
viscous dissipation. Moving outwards, a point on the disk 
will be able to see some of the star, although it is still par- 
tially blocked. The temperature on the concave side far 
from the warp will be between that of a flat passive disk 
which experiences no shadowing and a fully shadowed 
disk which is only heated by viscous dissipation. As the 
warp grows the temperature on the convex side increases 
and on the concave side it decreases. In limit of a very 
large warp the concave side will be completely shadowed 
while the convex side will be much warmer than a flat 
passive disk. 

Figure [5] compares the SED of the inner warp disk to 
that of a flat passive disk and a disk heated by viscous 
dissipation. The warped disk signiflcantly departs from 
the flat passive disk at short wavelengths because of the 
high temperature on the convex side while its SED is 
lower at long wavelengths because of the shadowing on 
the concave side. The warp drops below the midplane 
for 90 < 6* < 270° (the left side of the top half of the disk 
shown in figureflj)) and only half of the disk is shadowed. 
Since at most half of the disk can be shadowed, the flux 
always lies above a flat accretion disk even for very large 
warps. 

Figure |6] shows the variation in the SED as the warp 
moves around the star. The short wavelength flux 



changes as the convex side rotates more directly into 
the hne of sight while the long-wavelength flux does not 
change because this part of the disk is essentially flat and 
the SED of a flat disk does not depend on azimuthal view- 
ing angle. If the warp does not occult the star (left panel 
of figure [6|) there is a small change in the flux. View- 
ing the disk at a high enough inclination (right panel of 
figure |6]) changes the form of the SED variations. The 
infrared variations are smaller when the disk is viewed 
edge on than when it is seen closer to face on. This 
is because when the warp on the near side of the disk 
blocks the star, the hot side of the warp on the far side 
of the star is directly visible. These two effects work in 
opposite direction; blocking stellar flux reduces infrared 
flux, while seeing the convex side of the warp increases 
the flux. This can cause the change in optical flux and 
infrared flux to be in opposite directions, while for the 
more face-on case the flux variations are in the same di- 
rection at all wavelengths. 

The type of warp described here can be created by a 
number of different physical mechanisms. One possibility 
is a (sub) stellar companion embedded in the disk. As 
with the middle warp, if the companion is not coplanar 
with the disk then the gravitational perturbation from 
the companion can drive material out of the midplane. 
If the companion is located close to the star then this 
disturban ce will be at the inner edge of the disk. In the 
models of iLarwood fc Papaloizoul ()1997l ). for small mass 
ratios (q=0. 01-0.1), the warp changes the scale height 
of the disk by A/i/r w 0.1, similar to the scale height 
considered here. 

For the T Tauri star AA Tau, a warped disk caused 
by a tilted stellar magn etic field is beheved to explain 
the observed variations (jBouvier et al.l 120031 ). The stel- 
lar magnetic field is misaligned with the disk and the fiow 
of material onto the magnetic field warps the disk. This 
warp would occur near the corotation radius where the 
magnetic field from the star truncates the disk, where 
the dynamical timescale of the disk will only be a few 
days. Interferometric observations generally find that the 
hottest circumstellar dust is at distances consistent with 
the dust sublimation radius, much la rger than the coro- 
tatio n radius for higher mass stars (|Milan-Gabet et al.l 
|2007() . The exact location of the dust sublimation radius 
depends on the size of the dust grains and the geometry 
of the inner disk. If the disk is close to fiat, as is assumed 
in our models, or consists of large grains, then it may be 
possible for the dust to extend close enough to the star 
to trace a warp caused by a tilted magnetic field. 

The timescale for the inner warp to rotate around the 
star depends on the process originally responsible for cre- 
ating the warp. If the motion is precession driven by a 
companion on a misaligned orbit then the period wi ll 
be 1.5 years for a warp at 0.1 AU (|Caprani et al.l [20061 ). 
A warp created by a tilted magnetic field will rotate at 
the same rate as the star. Typical rotation periods for 



pre-main secjuenc e stars are a few days to a week (e.g. 

ICohen ' 

ations. 
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ICohen et al1l200l . 



which would produce rapid flux vari- 



A variable stellar accretion rate could also change the 
structure of the inner disk. The height of the inner rim is 
set by the incident luminosity, w hich is a combination o f 
stellar and accretion luminosity (jMuzerolle et al.ll2003D . 
Increasing the accretion luminosity will raise the temper- 



ature of the inner rim causing its scale height to increase. 
If the accretion luminosity is comparable to the stellar 
luminosity then a rapidly varying accretion luminosity 
could lead to a rapid change in the inner rim height. The 
accretion rate onto young stars varies r apidly on short 
timescales (H artigan et al.l [1991 : GuUbring et al.l [19961. 
If the accretion luminosity was not evenly distributed 
around the star, but instead confined to a hot spot then 
this may lead to preferential heating and growth of part 
of the inner disk. Growing part of the inner disk creates 
a non-axisymmetric structure that may be comparable 
to the inner warp structure, producing similar changes 
in the SED. 

The growth of an inner warp can create an anti- 
correlation in the disk flux as it did with the middle 
warp. In flgure[7|we show a model for g — 0.005 — 0.1(/i = 
0.005 - Q.OIAU) at i = 60°. As the warp grows the short 
wavelength flux increases while the long wavelength flux 
decreases due to the shadowing effect discussed above for 
the middle warp. The large temperature change in the 
inner disk leads to the signiflcant variation in the near 
and mid-infrared flux. This flux variation would be easily 
observable in the near and mid-infrared. 

The timescale to change the structure of the disk is the 
dynamical timescale, equal to the keplerian period, which 
is roughly one week at 0.1 AU. After the dust moves 
within the disk, its temperature will readjust to its new 
orientation relative to the radiation held on a thermal 
timescale. For the surface layers of the disk responsi- 
ble for the mid-infrared emission this timescale is a few 
seconds, effectivel y instantaneous compared to the dy- 
namical timescale (iChiang fc Goldreichlll997[ ). We com- 
pute warps of different sizes independently of each other, 
which is the same as assuming that as the warp grows 
its temperature responds instantaneously to the chang- 
ing intercepted radiation held. These rapid variations in 
the scale height could be caused by the motion of a com- 
panion misaligned with the disk, as in the middle warp 
([Fragner fc NelsonI [20091) . 

2.5. Spiral Wave 

We can also use the framework of TB96 to derive the 
SED of a simple spiral wave. Spiral waves are a common 
feature in the disks surroundin g stars with plan e ts or 
a ste llar companion (Gold rcich fc Tremaind Il979t iKlevI 
Il999| ). The shocks created by the passi ng density waves 
can change the scale height o f the disk (jEdgar fc QuillenI 
l2008t iBolev fc Durisenll2006D causing heating and shad- 
owing effects similar to a warp. We use a simple parame- 
terization of a spiral wave as a modiflcation of a flat disk 
to study its effects on the SED. The position of the spiral 
wave is given by: 



' min — sw 



(1 + nO) 



(11) 



The value of n sets the maximum radius of the wave 
and we keep n flzed at 0.8/(27r ) to match more detaile d 
descriptions of a spiral wave (jOgilvie fc Lubow[ l2002| ) . 
The parameter rmin-sw sets the minimum radius for the 
center of the wave, which is taken to be slightly larger 
than the innermost radius of the disk. Our wave only 
runs from < 9 < 360°, although it could continue to 
wrap around the star until it reaches the outer edge of 
the disk. The height above the midplane is given by: 



Hr,0)=g[l- — ]r. 



mtn — sw 



exp 



(12) 
The shape of the wave is a gaussian with width a,, cen- 
tered at rsw The function (l — -1^) modifies the height 
of the wave so that it decreases as it extends around the 
disk. Spiral waves grow weaker as they get further from 
their point of origin and the parameter m can be used 
to modify how quickly the wave height decreases. In our 
fiducial model the wave disappears at 27r (to=1). The 
factor g serves a similar purpose as before; it sets h/r for 
the wave. Although this parameterization allows us to 
modify many details of the spiral wave we generally keep 

771 ^ 1, (Tj. = U.lr„iin, "Tmin — sw ^ "^min +2*0"^, 77 := U.o/zTT 

and modify g. 

The temperature structure of spiral disks with g = 
0.05,0.1 is shown in Figure [H On the side of the wave 
facing the star the disk reaches a very high temperature, 
while directly behind the wave the temperature drops. 
Further out from the wave the disk is only partially shad- 
owed and the temperature lies between that of a flat pas- 
sive disk and the base flat accretion disk. In the models 
shown in figure [5] the radius of the wave increases with 
9 while the height decreases, hence at 6* = 180°, shown 
as the grey line, the temperature is lower on the front 
side of the wave and higher behind the wave than for 
9 = Q° . Comparing the two models, we find that as 
the size of the wave increases the temperature of the side 
of the disk facing the star increases while the outer disk 
becomes more shadowed. 

The resulting SED for g = 0.05 is shown in Figure [51 
It rises above the flat passive disk shortward of log A = 
1.0 while longward it drops significantly below the flat 
passive disk. The short wavelength increase over a flat 
passive disk is smaller for the spiral disk than for the 
inner warp because less of the disk is heated up to a high 
temperature (compare Figure [5] and [S]). The decrease 
at long wavelengths is much larger because in the outer 
disk every 9 is exposed to a wave inside its radius. For 
the warp only half of the disk is shadowed at one time, 
while the spiral wave can shadow all of the disk (figure [T]). 
Comparing the temperature of the outer convex side of 
the inner warp disk with the 9 = 180° spiral wave we 
can see that the temperature is lower for the spiral wave, 
which will result in a lower flux at longer wavelengths. 

An embedded massive p lanet within the dis k 
will create a spiral wave (jPapaloizou et al.l I2007D . 
These waves drive shocks which can substantia lly in- 
creas e the scale height of t he disk (Bolcv & DurisenI 
[2006'; 'Edgar & QuiUcn' [200 8[) . up to A/i/r « 0.6 
([Bolev fc Durisen .2006}) . The shadow cast by the 
planet on the outer disk could be imaged in scattered 
light, but would be difficult to observe in the infrared 
(| Jang-Condelll [20091) . The presence of massive planets 
have often been invoked to explain the lack of short- 
wavelength excess within young circumstellar disks, and 
spiral waves are a natural extension of this situation. A 
spiral wave created by a planet embedded in the disk will 
stay fixed relative to the planet and will travel around 
the star at the orbital frequency of the planet. For a 
planet within 0.1 AU of a star this can approach one 



week, which is much shorter than the precession of the 
warped disks. We can demonstrate the change in SED 
as the planet orbits around the star by observing the spi- 
ral wave at different azimuthal angles, shown in figure [HI 
The flux changes by < 0.05 dex and the change in the 
short and long wavelength flux is correlated, when the 
disk is viewed close to edge on. The small change in flux 
is due to the spiral wave being closer to an axisymmetric 
structure than either the middle warp or the inner warp. 
If the disk is highly inclined then the spiral wave will oc- 
cult the star and this occultation will dominate the vari- 
ations in the infrared (right panel of figure [S|). Here the 
infrared variations are substantial and strongly depen- 
dent on viewing angle, but are mainly due to changing 
the underlying flux from the star rather than changing 
the flux from the disk. 

Varying the wave height will change the shadowing 
of the disk and can produce an anti-correlation in flux 
f Figure [TUj) . as in the previous models. In these mod- 
els {g = 0.005 - 0.1, h = 0.0006 - 0.013 AU) we view 
the disk along a ~ 180° in order to directly observe 
the warmest part of the wave and produce the largest 
mid-infrared variations. As the wave grows the short- 
wavelength (A < 10/im) flux does not change as much as 
the inner warp because only a small part of the disk will 
be heated as it grows. The extra heated part of the spi- 
ral disk is confined to the inner edge of the wave, whose 
position does not change. When the inner warp grows, 
more of the disk is directly exposed to the star. This 
means that not only does the temperature at one point 
rise, but the area of the disk that departs from a flat disk 
is larger. The spiral disk only has the effect of increasing 
the temperature at a point, without increasing the area 
over which the disk is warmer than a flat disk. Beyond 
10/im the disk flux drops dramatically with an increase 
in the height of the wave. This is because all angles 9 
are subject to being shadowed; for the inner warp only 
—90° < 9 < 90° will be shadowed. This changes can 
reach up to ^ 0.2dex at 30/im. Changes in the height 
of the wave may be cause d by the nonlinear sho ck as 
it passes through the disk (jBolev fc Duriseiil 120061 or it 
may be caused by localized accretion hot spots heating 
different portions of the wave, as discussed for the warped 
disk models. 

Our spiral wave model can also be used to simulate 
other disturbances in a disk that appear as a wave. One 
such situation is a thermal instability in which a verti- 
cal disturbance grows i n to a wave that slowly travels 
inward toward the star (jWatanabe fc LinI 120081 ). In ac- 
cretion disks the temperature and structure of the disk 
are closely connected since the scale height depends on 
the disk temperature at a given radius. If the disk is 
flared, it is more directly illuminated by the stellar radi- 
ation field and hence is warmer at a given radius than if 
it were flat. If a small disturbance raised a piece of the 
disk above the midplane it would heat up as it became 
more directly illuminated. The frontside of this distur- 
bance will grow because it is warmer, while behind the 
disturbance the disk is shadowed and it will shrink. As 
a result, t he disturbance wil l propa gate towards the star 
as a wave. iWatanabe fc LirJ (120081 ) find that this type of 
instability can develop in the outer disk and is long lived 
enough to travel within a few AU of the star. 

We treat the instability as a spiral wave whose radial 



position and height do not change as a function of 9 
within the disk (n = 0, m = 0), and compare three 
locations and heights of the instabihty (r=0.5,0.4,0.3 AU 
with g=0. 05, 0.02, 0.01 resp ectively) that are typ ical of 
thermal instability models (jD'Alessio et al.1 11999"). The 
results, shown in figure [TT] have a pivot point in the SED 
at much longer wavelengths than the other models. At 
>3/im the flux change is ~0.1 dex until the far-infrared 
where the flux change can be much larger. In the models 
with g = 0.05, 0.02 the entire outer disk is shadowed and 
there is very little change in the far-infrared flux between 
these models. The flux from the outer disk only changes 
for the g=0.01 model where less of the disk is in shadow, 
which increases its flux. These waves travel inward on a 
thermal timescale, which is on the order of years. 

3. COMPARISON TO OBSERVATIONS 

3.1. LRLL 31 

iMuzerolle et al.l ()2009[ ) suggested that the mid-infrared 
variability of LRLL 31 may be due to a non-axisymmetric 
structure in the inner disk, and with our simple 
models we can more directly test this hypothesis. 
The observation s of L RLL 31, described in detail by 
IMuzerolle et al.l (|2009f) and summarized h ere, include 
multiple epoc hs with the Spit zer IRAC (iFazio et al 



2003), MIPS (|Riekeet al.ll200l and IRS (IHouck et al 



2004) instruments. Originally identified as variable based 
on two earlier IRAC and MIPS surveys of the field, it 
was followed up with further imaging and spectroscopic 
observations in the fall of 2007 and the spring of 2008. 
The entire set of data are shown in SED form in Fig- 
ure [121 IR-S spectra taken a week apart (October 9,16 
2007) show that the 5-8/im flux decreases by 0.2 dex, 
while the 12-40/im flux increases by a similar amount. 
After four months (February 24, 2008) the 5-8^m flux 
has risen, with a corresponding decrease in 12-40/im flux, 
close to its original state. Another IRS spectra taken 
a week later (March 2, 2008) show the disk still in its 
high 5-8/im state with the same anticorrelated change 
in the short and long-wavelength flux. IRAC observa- 
tions taken from the G TO (Lada ct ah 2006 ) and C2D 
p0rgensen. et al.ll2006t ) surveys (February 11, 2004 and 
September 8, 2004 respectively) also show substantial 
changes consistent with the IRS spectra. The MIPS 
24/im measurements include observations from the orig- 
inal GTO and C2D (RebuU et al. 2007) maps (Febru- 
ary 21, 2004 and September 19, 2004 respectively), as 
well as more intensive observations from our GO pro- 
gram (PID 40372). The MIPS foUowup includes five 
consecutive days (September, 23-27, 2007) and two fur- 
ther maps the following spring (March 12,19, 2008). The 
24/im fluxes change by a few tenths of a magnitude on 
weekly to monthly timescales. The uncertainty in the in- 
frared observations is < 5% while the changes range from 
20-60% suggesting that the variations being observed are 
real. In Figure [12] the October 9, 2007 spectra is shown 
with error bars to demonstrate the typical uncertainty 
in the spectra, which is smaller than the observed flux 
changes. A more detailed description of the data reduc- 
ti on, as well as the uncert ainties in the data, is presented 
in IMuzerolle et al.l ()2009l ) . All of our observations show 
variations at all timescales over the entire observed wave- 
length range. 



The low 5-8/im flux seen in the October 16, 2007 obser- 
vation is difficult to explain with the standard axisym- 
metric puffed inner rim model. In that case, the amount 
of dust emission is set by the scale height, which in turn 
depends on the stellar plus accretion luminosity. In the 
limit of low Lace and/or high _L», as is the case for LRLL 
31, the irradiation of the inner rim cannot change signifi- 
cantly, and its emission should always be considerably in 
excess of the photospheric level. Disk models with large 
scale asymmetries such as warps or spiral waves, coupled 
with highly inclined viewing angles, offer a way around 
the problem. 

We take LRLL 31 to be a G6 star, based on NIR spec- 
tra (jMuzeroUe et al.l [20091) . with a ra dius of 2.9i ? and 
a mass of I.8M0 derived using the iSiess et aLJ |2000[ ) 
isochrones. It has high polarization and its extinction 
is roughly Ay = 10. When presenting the observa- 
tions belo w, all o f the data have been dereddened by 
Ay = 10. iDahml ()2008D find the accretion rate in the 
range 10~^ — lO~^M0j/r~^ based on optical veiling and 
various emission line diagnostics. The inclination of the 
disk has not been directly measured, but given the high 
extinction and high polarization, which are signs that 
some of the light passes through the disk, we focus on 
inclinations close to edge on. 

Given the modeling results presented earlier, we focus 
on changing the height of the warp/wave. The middle 
warp, inner warp and spiral wave all produced a wave- 
length dependence similar to LRLL 31, where the short- 
wavelength flux would increase while the long- wavelength 
flux would decrease due to the shadowing of the outer 
disk. Figures [131 (HI and [151 show the middle warp, in- 
ner warp and spiral wave close to edge on with varying 
heights along with the observations of LRLL 31. An 
inclination of i=85° was chosen for these three models, 
while different values of a were chosen to maximize the 
flux change. The middle warp displays a similar wave- 
length dependence as the observations, but is unable to 
match the size of the variations at 3.6,4.5/im. Also the 
location of the pivot in the SEDs is at a longer wave- 
length than observed. The inner warp does a better job 
of fitting the size of the variation as well as the slope of 
the IRAC and short-wavelength IRS data. For the inner 
warp we consider two inclinations that are both close to 
edge-on but produce different SEDs. The difference in 
the two edge-on inclinations comes in the optical where 
the i=85° disk occults the star while the i=95° disk does 
not. This model also shows a pivot point in the SED, but 
it is also at a longer wavelength than observed. For the 
spiral wave model the pivot point moves to a wavelength 
that is close to 8/im but the short wavelength flux vari- 
ation is too small to match the observations, in particu- 
lar the 3.6,4.5/im photometry and the spectra with the 
lowest 5-8/im flux. The timescale for the warp to change 
height, which is the dynamical timescale, approaches one 
week at 0.1 AU around LRLL 31. This is short enough 
to match the rapid variations seen in the mid-infrared. 
We do not consider the thermal instability because the 
timescale is much longer than one week. 

Based on our simple models, we find that a disk with 
an inner warp whose height rapidly varies is able to best 
explain the observations of LRLL 31. It can fit the flux 
change at A < 8/im, it includes a pivot point in the SED, 



and the timescale is consistent with the observations, but 
the simpUcity of the model prevents it from fitting in 
detail. More detailed radiative transfer models would 
help to further constrain the shape of the perturbation. 
Future work should also tie these parametric models to 
specific physical mechanisms in order to provide predic- 
tions for the detailed time-dependent behavior of the flux 
changes. 

3.2. Application to Other Variables 

In the previous models we considered a star whose mass 
was on the high end for a typical T Tauri star and a lower 
mass star with a lower effective temperature will not heat 
the disk to as high of a temperature at the same dis- 
tance from the star, which will result in a weaker excess 
for a given value of h/r. As the disk precesses, only the 
orientation of the disk relative to the observer changes, 
which does not depend on the central star or the tem- 
perature at a particular location within the disk. The 
change in infrared excess with precession will then be 
similar over a wide range of masses. The same is true 
if we vary the scale height of the disk. In figure [12] 
we show an inner warp seen at i = 60° around LRLL 
31, a typical T Tauri star (Te//=4000 K, M^^IMq, 

i?, = 2Rq,M = 5 * lO-^Mgyr-i, R^,n = Si?*) and a 
brown dwarf (Te//=3000 K, M,=O.O8M0, R^=0.7Rq, 
M = lO^^Afgyr^^, Rmin ~ 8i?,) with a height varying 
from g = 0.01 — 0.1. The accretion rates were chosen 
to decrease with mass, but they only contribute signif- 
icantly to the heating of the outer disk where our as- 
sumption of a flat disk severely underestimates the flux 
from typical disks. The change in the mid-infrared flux 
around all three stars is almost identical. While the exact 
structure and temperature of these three disks differs, as 
the mass decreases the temperature of the disk decreases 
and the physical distance between the disk and the star 
is smaller, the change in the warp is the same resulting 
in similar variations in the infrared flux. Changing from 
an inner warp to a middle warp or a spiral wave will also 
show a similar lack of dependence on the central source 
as the height of the warp/wave varies. If the structural 
change in the inner disk does not depend on the mass of 
the central star then such a change in the disk structure 
would be easily detectable in the near and mid infrared 
flux from many young stellar objects. 

Since changing the spectral type of the central star 
does not have a significant effect on the infrared variabil- 
ity presented here, these models can provide insight into 
other variable objects. iMorales-Calderon et al I (|2OO90 
use IRAC to survey 3-8/im variability among young stel- 
lar objects within IC 1396A. They find that variabil- 
ity occurs in more than half of the young stellar ob- 
jects, typically showing colorless variations with am- 
plitudes of 0.05-0.2 magnitudes on timescales of a few 
days to a week. This is consistent with our models 
of a growing inner warp, which when the warp does 
not occult the star, can produce large variations in flux 
without a large change in the c olor in this wavelength 
range. iHutchinson et al.l (|1994D find two stars whose 
BVRI fluxes and 4.5,10/xm fluxes vary in opposite di- 
rections while the near- infrared fluxes stay constant. 
iJuhasz et al.l ()2007D find similar results for SV Cep where 
the V band and 3.6/xm flux are anti-correlated. They also 



find that the V band and lOOyitm flux are correlated with 
each other, while there is very little change in the flux 
at 25/xm. This is qualitatively similar to our model of a 
growing inner warp seen ati = 85° f Figure [T4|) . As the 
warp grows the star becomes blocked, dropping the op- 
tical flux, and more of the outer disk is shadowed, drop- 
ping the long wavelength flux, producing a correlation 
between the optical and lOO/im flux. Also the larger warp 
emits more flux in the mid-infrared, producing an anti- 
correlation between the optical and 3-10/im flux, while 
the near-infrared bands stay roughly constant. 

UX Ori stars are a type of optical variable that exhibit 
occasional dimming and reddening of their optical colors 
due to occultations of the star by clumps of gas and dust. 
These clumps are thought to arise from the inner edge 
of the circumstellar disk, which can oc cult the star since 
the sy stem is seen very close to edge on (jDuUemond et all 
l2003f l . In some models of the infrared excess from these 
stars it has been assumed that the disturbance at the 
inner edge does not cover enough area to substantially 
shadow the disk ([Pontoppidan et a"Lll2007[ ). Our models 
show that there can be significant mid-infrared variations 
even if the perturbations at the inner edge of the disk do 
not shadow the entire disk. 

Ongoing studies of mid-infrared variability will help to 
determine how common this type of variability is among 
young stellar objects. Multiple programs in progress 
gathering multi-epoch Spitzer data at 3.6 and 4.5/im will 
constrain the frequency of mid-infrared variability among 
T Tauri stars. Given the wavelength coverage of these 
variability studies, and the expected mid-infrared vari- 
ability due to non-axisymmetric variations in the inner 
disk, any changes in the structure of the inner disk will be 
readily observable. Using these observations to look for 
periodicity, which is a sign of a perturbing companion, 
or a correlation with accretion rate will help to constrain 
the physical processes that influence the inner structure 
of circumstellar disks. 

4. CONCLUSION 

We have used simple geometric models to investigate 
the effects of non-axisymmetric structure in a circumstel- 
lar disk on the system's SED. By simplifying the models 
we have been able to investigate a wide range of geome- 
tries including a warp in the middle of the disk, a warp 
at the inner edge of the disk and a spiral wave. These 
models show an anticorrelation between the flux at short 
and long wavelengths because of variable shadowing of 
the outer disk, in qualitative agreement with many ob- 
servations. However, precession or corotation of these 
structures produces a negligible change in the infrared 
excess amplitude. Only by changing the scale height of 
the perturbations, as may be produced by the dynamical 
interactions or variable illumination of disk structures, do 
we find significant flux variations as observed. Further 
observations spanning a range of timescales and covering 
multiple objects will help to discriminate between models 
as well as indicate the frequency of such structural varia- 
tions. The similarity in the size of the infrared variations 
across a wide range of masses also indicates that the type 
of non-axisymmetric variability explored here would be 
detectable around many young stellar objects. Numer- 
ous observations of mid-infrared variability around young 
stellar objects are consistent with our models. 
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Fig. 1. — Schematic drawing of our three simple geometric models. 
(c) the spiral wave all viewed edge on along i = 90°, a = —90°. 
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The three models are (a) the middle warp, (b) the inner warp and 
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Fig. 2. — Middle Warp disk. In the left and middle panels are the temperature structures of a disk with a warp that reaches its maximum 
h/r of g = 0.05,0.01, respectively, at r = 0.5AU with a flat disk extending beyond the warp to 100 AU. The dark solid line is the side of 
the disk facing the star, while the grey solid line is the side of the warp facing away from the star. On the right is the SED for a disk with 



an middle warp that grows to <; = 0.05 at Ku_ 



0.5 AU and a flat extension to the disk that extends to 100 AU. The disk is viewed 



along i = 60° , a = 0° . Also shown are a flat accretion disk, a flat passive disk and a blackbody at T^f f = 5700i<' for comparison. 
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Fig. 3. — SED of a precessing disk with an middle warp extending to Rwarp = 0.5 AU with g = 0.05 and a flat extension seen at an 
inclination of 60° (left panel) and 85° (right panel). From dark to light lines the disk rotates through 180°. The inset shows a closeup of 
the models near lO^m. The dotted line is the photosphere. The change in the infrared flux is small, but it increases as the disk is viewed 
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Fig. 5. — Inner warp disk. The left and middle panels show the temperature distribution for a disk with an inner warp of g = 0.05 and 
g = 0.01 respectively. The dashed line is the temperature of a disk heated solely by viscous dissipation (M = 10~^MQyr^^), the dotted 
line is a perfectly flat passive disk, while the solid (dark and light) lines are the warped disk (convex and concave sides respectively). The 
maximum temperature of the warp, and the shadowing of the outer disk changes with the size of the warp. On the right is the SED 
for a warped inner disk with g = 0.05 seen at i = 60°, a = 0°. For comparison a flat passive disk, a disk heated by viscous accretion 
(M = 10~^M(7)yr~^) and the photosphere are also included. 
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Fig. 10. — SEDS for a growing spiral wave. The disk is viewed at i = 60°, a = 180° and the models vary g from 0.005 to 0.1. The change 
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Fig. 12. — Observations of LRLL 31. Black circles are the photometry from 2MASS, two epochs of IRAC photometry and 10 epochs of 
MIPS 24/xm photometry. Colored lines are the four IRS spectra of LRLL 31. The dark blue and light blue were taken in October 2007 and 
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representing the stellar photosphere for comparison, and the observations have been scaled to match the photospheric luminosity at J band, 
assuming the J band flux is completely photospheric. 
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Fig. 13. — SEDs for a growing middle warp (g = 0.01 — 0.1 from dark to light lines) viewed at i = 85°, a = 180°. On the left is the entire 
model, while on the right is a close-up of the model with the observations. The observations have been scaled to match the J band flux of 
the model, assuming the J band observation is completely photospheric. 
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Fig. 14. — SEDs of LRLL 31 with a growing inner warp. The models vary g from 0.005 to 0.1 (from dark to light lines). The models on 
the left are seen at an inclination of j = 85° , while on the right are the models seen at i = 95° . The top panels show the entire models 
while the bottom panels zoom in and include the observations. Only models for g = 0.005,0.02,0.04,0.06,0.08,0.1 are shown for clarity. 
The observations have been scaled to match the J band flux of the model, assuming the J band observation is completely photospheric. 
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Fig. 15. — SEDs for a growing spiral wave around LRLL 31. The disk is viewed at i = 85°, a = 180°. The models vary g from 0.005 to 
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Fig. 16. — SEDs for an inner warp surrounding LRLL 31, a typical T Tauri star and a brown dwarf seen at i = 60°. For each star the 
models vary g from 0.01 to 0.1(dark to light grey solid lines) and the stellar photosphere is shown with the dotted line. The change in 
mid-infrared flux docs not depend strongly on the central star. 
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APPENDIX 

5. OUTER WARP + FLAT EXTENSION 
5.1. Adding a Dependence to Outer Warp 

TB96 describe how to calculate the SED for an outer warp seen at various inclination angles. Their general method 
for calculating the SED also includes a dependence on azimuthal viewing angle, although their detailed treatment of 
various occultation effects (the star blocking the far side of the disk, the disk blocking the star, etc.) does not include 
this dependence. Since these disks are non-axissymetric the SED can depend substantially on the azimuthal viewing 
angle, a, of the observer. In this section we describe how we have added this dependence into the equations of TB96. 
We do not include a detailed description of the derivation of the equations, but merely state most of them and some 
of the geometric logic behind their modification. 

The first modification is to equation (6) of TB96, which describes the calculation of the flux from the disk based on 
the temperature of the disk. This equation assumes that the disk is viewed along a = so that the upper disk from 
to 7r/2 looks the same as the upper disk from 37r/2 to 27r. When a ^ this symmetry is broken and the individual 
components of the disk must be considered. While the temperature of a concave, or convex, piece of the disk does 
not change with viewing angle, the orientation of each of the concave, or convex pieces, changes and must be treated 
separately. For the outer warp this becomes: 



px/2 , s i-it/2 



^u.u — J^ Jq ^[-^ concave) Jup i Jq ^ \-^ convex ) .flow i J7r/2 \ co7icave)jlow i J7r/2 \ co7ivex ) Ji 



rRdisk 



l-3Tr/2 „,„ N r c37r/2 



Jr J IT ^\-^ convex ) J up i J^ J^\ J- concave) Jlow i jZ-k/2 \ concave) Jup i jZ-k/2 \ convex) Jlow 1*^7 

The next change comes in the appendix to the functional form of the parameter C. The function C is used to deflne 
the points that are along the line of sight with the star. If the line intersects the star then we need to worry about 
whether the disk blocks the star or the star blocks the disk. If this line does not intersect the star then the disk cannot 
block the star and the star cannot block the disk. The definition of C changes from C — rsin6' to C = rsin(6' — a). 

Also the radial part of the deformation used in equation (A6) of TB96 is taken to be H{r) = gRdisk ( -^. — ) cos a. 

For a ~ 7r/2 the disk along the line of sight is flat and the radial part of the height will remain at zero, while along 
a — TT the disk curves below the midplane as expected. 

The flnal change comes when calculating the stellar flux. In equation (A12) of TB96 we take h{r,9) to be h{r,a) 
since this represents the part of the disk that will block the star. As the azimuthal angle increases the disk blocks less 
of the star because the height of the disk is smaller. We make more changes to how the stellar flux is calculated, which 
are described below, but when it comes to the occultation of the star by the disk this is the only change. 

In the end we are able to run our models from < i < tt and — 7r/2 < a < 7r/2. Due to the symmetry of the disk 
this covers all possible viewing angles allowing us to accurately model the precession of the warp, as well as observe 
the warp from an arbitrary angle. 

5.2. Flat extension of Outer Warp 

We have taken the outer warp model and added a flat extension beyond it in order to treat disks where the warp 
is not at the outer edge of the disk. The warp will shadow the outer disk, changing its temperature structure. For 
simplicity we assume that the outer disk is a flat blackbody. The temperature can be derived using the same formula 
as with the warped disk (Equation 6), but with different deflnitions for the integration boundaries. Half of the flat 
extension will be shadowed while half will not be shadowed. For the side of the flat extension that is beyond the part 
of the warp the goes below the midplane there is no shadowing of the disk and the integration ranges over: 

^min U 
^niax TT/Z 

(14) 

^min U 

Smax = arcsin(i?*/(i) 
d^ = r^ + h^ 

(15) 

For the part of the flat extension that lies behind the warp that stretches above the midplane the definition of 6min 
changes. In this case dmin is set by the angle between the warp and the point P(r, 9) in the disk. 

9RwarpCOs{9) \ 

(r - H ) ^ ' 

This takes into account shadowing of the flat extension to the disk due to the warp. Once the temperature structure 
has been determined the flux can be derived using the equation for the flux from a disk (Equation 5). 
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Fig. 17. — Schematic diagram of the limits on 5 for the convex side of the disk. 




Fig. 18. — Schematic diagram of the limits on 5 for the convex side of the disk. In this case the star is blocked by some of the disk. The 
shallow slope of the disk prevents all of the star from being seen. 

6. INNER WARP 

6.1. Temperature Profile of the Inner Warp 

In this section we describe the method for calculating the SED of a disk with an inner warp. In the text we laid 
out the basic equations from TB96 that are needed to calculate the temperature structure. As mentioned in the text 
the essential difference between the inner warp and outer warp comes in calculating Smax,Smin,£min,£max for each 
point P{r, 9), which are used in equation 6. From here the disk is split into two sides that are treated separately. The 
convex side is the side of the disk that faces the star on the inner edge and receives the most direct heating from the 
star. For this side, the integration ranges over: 



-rmn 



==0 

^max — IT I ^ 

(17) 
5min = — arctan(9/i/3r) 
^raax = arcsin(i?,/d) 
d? ^r"^ + h^ 

(18) 

Figure [T7] demonstrates the limits on 5 for the convex side of the disk. The definition of 5min , demonstrated in 
figure [TU comes from the inner disk blocking light from the top of the star. The inner disk will limit the field of view 
of the point P{r, 6) as it looks toward the star. Traveling out from the star, less of the star will be seen by the disk 
because the shallower slope of the disk will cause more of the star to be blocked. In the limit of a flat disk far from the 
star 6min approaches zero and the disk can only see half of the star. If the point P{r, 9) on the disk is close enough to 
the star then the disk can see all of the star and 5„iin — ~Smax- The limits on e assume that the scale height of the 
disk does not change across the face of the disk, which will be an accurate approximation far from the star. 

The concave side of the disk is the side that does not directly face the star. Since it does not face the star much of 
the inner parts of the disk will be blocked by the warp and will only be heated by viscous dissipation. The condition 
for the point P{r, 9) on the concave part of the disk to see any of the star is: 

If the point P{r, 9) meets this condition then this point can see some of the star and Smin becomes f fig [TO]) 

5min = arctan( — -^^^^ — ) (20) 

The rest of the limits stay the same. In the limit of a perfectly fiat disk Smin approaches and the point P{r, 9) 
is irradiated by only half of the star. For a large warp the only heating by this side of the disk will be from viscous 
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Fig. 19. — Schematic diagram of tlic limits on 5 for the concave side of the disk. 



dissipation because the warp will block the star over most of the disk. 



6.2. Calculating the SED for an Inner Warp 

In this section we describe the procedure for converting the temperature structure into a spectral energy distribution 
(SED). From TB96, the flux emitted by the disk is 



i^..u= / / B,(T(r, 0))d5nd-u (21) 

J J disksur face 

In this case u is the vector along the line of sight to the observer from the center of the star and nd is the normal 
to the disk at the point P(r, 9). The vector u can be defined in terms of the azimuthal and polar angles to the line of 
sight, a and i respectively. The area of the disk along the line of sight is given by 

dSriA • u = r [(iflsine - |^cos6') cos a sin z - (i||cos6' + ^sitlO) sin a sin z + cos i] drd9 (22) 

The angle a ranges from — 7r/2 to 7r/2 while the inclination i ranges from to tt. This covers all possible viewing 
angles of the disk, since the symmetry of the disk makes some viewing angles redundant. 

Splitting up the equation for the flux from the disk helps to make the problem simpler to understand and more 
tractable. It also fits with the fact that we do not need to calculate the temperature structure of the entire disk. The 
symmetry of the disk allows us the calculate the temperature of the convex and concave side from Q < 9 < tt 12 and 
then apply this temperature profllc to the rest of the disk. The integral is split into eight parts: 



F„ 



jj. jQ ^\-^ concave) J up > Jq ^[-^ convex ) flow > J ^ /2 -^\-^ concave) Jlow i J^ /2 -'^[-^ convex ) Jup 



cRdi. 



r3TT/2 j-.,rTi \ r r37r/2 „ .„ ^ ,. f2TT 



~^Jr J 7T ^\-^ convex ) Jup I j ^ ^\-^ concave) J low ~^ jr>^^ /2 ^\-^ concave) J up y J^^^ 12 ^\-^ convex) Slow V '^J 

where /(r, 9) = dSn^ ■ up(r, 9), and p{r, 9) is a binary function used to determine if the point P(r, 9) is visible to the 
observer. The integration is done over both the upper and lower sides of the disk in order to account for inclination 
angles greater than 90° where the lower half of the disk is visible. If the inclination is then the observer is face on 
to the upper half of the disk, which has both a concave and convex side. If the inclination is 180° then the observer 
is face on to the lower half of the disk, which includes both a concave and a convex side. Treating each quarter of 
the disk separately allows us to use the symmetry of the temperature proflle but still treat general azimuthal viewing 
angles. 

6.2.1. Calculating the value of p 

The above description sets out the basics for how to calculate the temperature structure and SED for a warped 
inner disk. Most of this is derived from TB96, which treated these situations generally enough to apply to any type 
of warp. The main differences between this inner warp and the outer warp from TB96 comes from the calculation of 
p{r, 9). This section describes the conditions used to calculate p(r, 9) for the particular warp used here. 

The flrst condition is that the observer is facing the point P{r,9). For inclinations less than 90° the observer will 
see mostly the upper half of the disk, while at inclinations greater than 90° the observer will see mostly the lower half 
of the disk. There are select inclinations close to edge on where at inclinations less than 90° some of the lower disk 
can be seen. For example, if figure [TH] had an observer in the upper left viewing the disk close to edge on they would 
be able to see some of the lower convex side that is illustrated in the figure. In general it can be determined if the 
observer is facing the point P{r, 9) based on the dot product nd • u. The normal is defined as extending on the upper 
side of the disk and the dot product will be greater than zero if nd and u lie along the same direction. Therefore, the 
upper part of the disk can be seen if the dot product is greater than while the lower part of the disk can be seen 
when the dot product is negative. 

Now we determine if the point P(r, 9) is blocked by either the star or the disk. First we consider whether the star 
blocks the far side of the disk. This applies for inclinations less than 90° where the star may block part of the upper 
convex side, as is demonstrated in figure [20l The limit at which this condition applies is given by 
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Fig. 20. — Schematic diagram demonstrating when the star starts to block some of the disk. 



TABLE 1 
Limiting inclination for occultation of disk 

BY STAR 



2R, 



5i?* 



lOiJ, 



0.005 


59.7 


78.2 


83.9 


0.01 


59.4 


77.9 


83.7 


0.03 


59.3 


76.7 


82.5 


0.05 


57.2 


75.6 


81.4 


0.07 


56.1 


74.5 


80.3 


0.1 


54.5 


72.8 


78.6 


0.3 


44.7 


62.3 


67.8 


0.5 


36.7 


53.1 


58.3 




Fig. 21. — Schematic diagram demonstrating the part of the disk that gets blocked by the star. 



tailZ;j„j — 



i?* COS i 



lini 



grmin COS a + R^ sin in 



(24) 



If the inclination is greater than ium then part of the disk is blocked by the star. In this case we can use the discussion 
of TB96 to determine if the point P{r,9) is blocked by the star (figure [H] and table [T]). Defining C = rsin(6' — a) 
the only time the star can block the disk is when C < R^ otherwise p{r,9) = 1. If C < i?* then p{r,9) ~ 1 only if 
rcosO < td cos 9d, where 9d = tt — arcsin(C/r£)) and ro is the positive root of the following equation: 
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Fig. 22. — Schematic diagram demonstrating the part of the disk that gets blocked by the star. 




Fig. 23. — Schematic diagram demonstrating when the disk is steep enough to block itself. 
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zm smi + 




C_ 



+ xm cos z = 



(25) 



where H[r) — grmin ( ''. ) cos a and M is the point on the northern hemisphere of the star in the plane (P,u,z) 
such that u is tangent to the star at this point: 



Xm = —\/Rl - C^ cos I 

TJM = C 

ZM = VR* - C^ sin i 



(26) 
(27) 
(28) 



Now consider inchnations greater than 90°. In this case the disk can stih be blocked by the star if the warp is small 
enough and the inclination angle is close enough to 90 degrees (figure [22| . This affects both the upper and lower part 
of the disk, but only from 9 E [a + 7r/2, a + 37r/2]. The condition for the point P to be hidden by the star is: 



\h{r,0)\< 



R^ cos(z — 7r/2) tan(7r — i) — r 
tan{'ii — i) 



(29) 



As above this only applies when C — Rt, sin(6' — a) < i?*. 

Another possibility to consider when the inclination is greater than 90° is that the warp on the lower concave side 
is steep enough that is blocks part of this side of the disk (figure B^ . This condition only applies to parts of the lower 
concave side, from 9 E [a + tt/2, a + 37r/2]. The condition is that 71 < 72 where 71 = i — 7r/2 and 



72 



h,n - h{r, 0) 
PM 



(30) 



In this case M is the highest point on the disk on the line of sight to point P. For r sin(7r — 0) > rmin the disk does 



block itself, but for r sin(7r — 0) < Tmin we have r„ 
distance between the point P and the point M (\/(7i^ 



and 0m = arcsin(r sin(0 — Q;)/r„). The quantity PM is the 



hY + {rm-rY). 
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Once all of the conditions have been considered the flux from the disk can be calculated using equation [531 These 
conditions would apply to any type of warp whose maximum height above the midplane occurs at the inner edge of 
the disk, as opposed to at the outer edge of the disk, regardless of the exact functional form of the warp (ie. power 
law vs. exponential). 

7. SPIRAL WAVE 

7.1. Temperature Profile of Spiral Wave 

The third type of disk that we attempt to model contains a spiral wave. As with the warped disks we follow the 
derivation of TB96 to derive the temperature structure and SED for this disk. The derivation for the temperature 
structure is very similar to that of the inner warp, only with slightly different definitions of the boundaries. For the 
part of the disk inside the wave, the disk is not blocked by the wave but the amount of the star seen can change. 
For points far from the wave, the disk is like a flat disk and Smm = 0. For points on the wave, as it rises above the 
midplane, more of the lower half of the star will become visible. How much of the lower half of the star is visible 
depends on the location and height of the point on the wave. In this case the lower limit on i5 is: 

Smtn = - arctan(/i/(r - r™i„)) (31) 

This limit will continue to increase until the point on the wave can see the entire star and then Smin = ^Smax- The 
other limits stay the same as in the previous models: 

^max ^ TT/Z 

(32) 
Smax = arcsin(i?*/(i) 
(f = r^ + h^ 

(33) 

For the parts of the disk behind the wave, some of the star may be obscured. In this case: 

Smin = arctan(/ism/(r - r^u.)) (34) 

hsw = grmin-swC^ " m9 /2Tr) 

"^ sw — '^min — swy^ ~r Tiu ) 

(35) 

This is similar to the concave side of the inner warp, where the warp can obscure part of the star. The only difference 
is that the maximum height of the disk does not occur at the inner edge of the disk, but instead occurs at the location 
of the spiral wave. When Smin > Smax then the entire star is blocked and that point on the disk is only heated by 
viscous dissipation. With these definitions and equation (6) we can calculate the temperature of the disk. 

7.2. Calculating SED of Spiral Wave 
The flux from the disk is given by: 

B{Tdisk)fup (36) 

There is no symmetry in the disk that allows us to split the disk into different parts, as with the concave and convex 
pieces of the inner warp. We also only consider inclinations less than 90°, where we only see the upper disk, since the 
lower disk will look the same as the upper disk. 

The one occultation effect we include is the blocking of the disk by the wave. For a — 7r/2 < < a + 7r/2, the near side 
of the disk, the wave can block the part of the disk that is at smaller radius than the wave. For a + 7r/2 < 9 < a + 37r/2, 
the side of the disk on the other side of the star from the observer, the outer disk can be blocked by the wave. This 
effect can become important for modest inclinations, given the typical wave heights we consider here. To determine if 
a point on the disk is blocked by the wave we first need to determine where the line of sight intersects the wave. This 
is illustrated in figure [24] and is given by: 

x = rsin^ = rmin-swC^ + n9M)sin{9M) (37) 

Here 9m is the azimuthal coordinate of the point where the wave intersects the line of sight. We assume that the 
point M lies at the peak of the spiral wave. This is only an approximation, although the narrowness of the wave make 
it an accurate one. The angle between the line connecting the points P and M and the midplane is 7 (Fig. [25]). When 
7 > tt/2 — i then point P{r, 9) is blocked. 
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Fig. 24. — Schematic diagram demonstrating the position of M and P in the disk when the wave can block part of the disk. 
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Fig. 25. — Schematic diagram demonstrating how the spiral wave can block part of the disk. 




Fig. 26. — Schematic diagram defining the angles ^,^ 



tan 7 



I'M 



ru - rp 



(38) 



We ignore occultation effects due to the star blocking the disk, which we did consider in the inner and outer warp 
model. Based on our experience with the warped disks and the typical dimensions of the disk, these are negligible 
effects that will only play a role very close to edge on. We also do not consider situations where the wave on the near 
side of the disk can block the far side of the disk. The exclusion of these two effects prevents us from considering the 
spiral wave at inclinations very close to 90°. 

8. STELLAR FLUX 

Next we consider the flux coming from the star. We follow a similar procedure as above where the equation for the 
stellar flux is modulated by a binary function (£(</>, '0)) which equals 1 when that part of the star is not blocked by the 
disk and it equals zero when the star is blocked by the disk. The angles 0, V'j shown in figure BBl are the azimuthal 
and polar angles of a point on the surface of the star relative to the center of the star and the z axis (the same z axis 
as for the disk). The x-axis of this coordinate system is in the same direction as the line of sight, and will differ from 
the X-axis of the disk by the angle a. The fiux from the star is: 



-F* — B,, 



,V)rfA 



(39) 



surface 



To determine the surface over which we integrate, we need to know the points of the star that are seen by the 
observer (ie. which side of the star is facing the observer). These points will be those that have u • dA > where 
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Fig. 27. — Schematic diagram showing the case where i < 90° and the disk can block the star. 



U : 



sinix + cosiz 



dA — R^ sin ipd^pdcjilcos (p sin ^x + sin (p cos ^y + cos tpz) 



(40) 
(41) 



The evaluation of e{(j), tp) will depend on the type of warp/wave and the orientation of the observer. First consider 
inclinations less than 90°. In this case the warp/wave streching above the midplane may block some of the star. The 
entire star will be blocked if the following condition is met: 



/i(r, a) — r tan(7r/2 — i) > i?^ 



(42) 



where r is location of the peak of the vertical disturbance and h is the maximum height of the warp or wave at the 
angle a. The exact value of r and h will depend on whether we are considering the outer warp, inner warp, or spiral 
wave (ie. r = Tmim h{r,a) = grminCOs{a) for the inner warp). This condition is illustrated in figure ?? for the inner 
warp, which is the disk that is the most likely to occult the star. 

None of the star will be blocked if the inclination is less than inm (discussed earlier). A generic version of the 
equation for ium that can be applied to all of the disks is: 



taniiim — 



i?, cos i 



lim 



h{r, a) + i?* sini;^ 



(43) 



When the inclination falls between these two limits only a fraction of the star is blocked. We can use the discussion 
of TB96 section A2.2 to determine if a point of the star's surface is hidden by the star 

The point Q is a point on the surface of the star that intersects the line of sight and the upper edge of the disk. If 
the point N{(j),ip) lies above Q then the observer can see this part of the star, otherwise it is hidden and e((/),V') = 0. 
The vertical coordinate of Q, zq, is the greatest root of the following equation: 



(1 + tan^ i)zQ ^ 2 tan 2 (/i tan z — r cos q;)zq + (/itanz — r cosa)^ — R^ + R^ sin^ ^cos^ (j) — 
If zn > zq then e{4>,ip) = 1 otherwise £(0, -0) = with zn being given by: 

zn = —R* sini/; sin cpsini + i?* cosi/" cosi 



(44) 



(45) 

For the inner and outer warp we consider inclinations greater than 90° where the disk can still block part of the 
star, although it is less likely because the disk curves away from the observer. This is illustrated in figure ?? for the 
inner warp, but can also apply to the outer warp in the limit that h goes to zero. The point at which the line of sight 
is perpendicular to the normal of the disk sets a limit to the distance z above the midplane that an observer can see. 
If this distance is less than the radius of the star then some of the star is blocked by the disk. The point D at which 
the the line of sight is perpendicular to the disk occurs when u • n = (figure ??). If this condition is met at a radius 
ru then a point N{(l), tp) on the stellar surface will be blocked if 



i?* cos ip > h{rB ,a) + tb tan(i — 7r/2) 



(46) 



If there is no point at which u • n = then ru — rmin and the same condition for being able to see the star is used. 
When this condition is met e{(j),ip) = 0, otherwise e{(p,ip) = 1. 

All of these different occultations are combined to determine e{(j), ip). The flux is determined by integrating over the 
entire surface and added to the fiux from the disk to create the observed SED. 



27 




Fig. 28. — Schematic diagram showing the case where i > 90 and the disk can block the star. 



